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We (re) consider in this paper the problem of tunnehng through an impurity in a quantum wire 
with arbitrary Luttinger interaction parameter. By combining the integrable approach developed 
in the case of quantum Hall edge states with the introduction of radiative boundary conditions to 
' describe the adiabatic coupling to the reservoirs, we are able to obtain the exact equilibrium and 

non equilibrium current. One of the most striking features observed is the appearance of negative 
differential conductances out of equilibrium in the strongly interacting regime g < .2. In spite of 
the various charging effects, a remarkable form of duality is still observed. 
, New results on the computation of transport properties in integrable impurity problems are gath- 

ered in appendices. In particular, we prove that the TBA results satisfy a remarkable relation, 
originally derived using the Keldysh formalism, between the order correction to the current out 
^ . of equilibrium and the second derivative of this current at T = with respect to the voltage. 
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I. INTRODUCTION. 



Tunneling experiments [|lj are one of the most efficient probes of the physics of Luttinger liquids, which is expected 
to describe the properties of one dimensional conductors. The case of spinless Luttinger liquids has already been 
' O ' extensively studied, both theoretically and experimentally, in the context of edge states in a fractional quantum Hall 
bar, where in particular, shot noise measurements have led to the observation of fractional charge carriers ||]. The full 
cross-over between the weak and strong backscattering regimes has also been studied ||-|^: it exhibits in particular 
a duality between Laughlin quasi particles and electrons that is the result of the strong interactions in the system, 
and, ultimately, of integrability. From a theoretical point of view, it must be stressed that crossovers in this type 
of problems can only be properly studied with non perturbative methods anyway. In fact, for the physics out of 
equilibrium, which plays a crucial role in the shot noise experiments for instance, numerical simulations don't even 
seem to be available. 

' Other one dimensional conductors where Luttinger liquid physics could be observed include carbon nanotubes 

or quantum wires in semiconductor heterostructures j7j . A key question for the latter examples is how to describe the 
\ application of an external voltage. In the fractional quantum Hall case, this turned out to be easy [|| because the left 
. and right moving excitations are physically separated (the Luttinger liquid is really the "sum" of two independent 
chiral ones), and put at a different chemical potential by the applied voltage. This will not be the case in a real 
quantum wire, where various charging effects have to be taken into account. 

The matter led to some active debating |^,^ , and now seems quite settled. We follow here the approach of , which 
, easily allows the inclusion of an impurity. We thus consider a gated quantum wire coupled adiabatically to 2D or 3D 
• reservoirs. As in Landauer-Buttiker's approach for non interacting electrons pO| , p^ , these reservoirs are assumed 
P5 ' to be "ideal" , and merely are there to inject bare densities of left and right movers in the wire. The interactions in 
O the wire lead to the appearance of a non trivial electrostatic potential, and, in turn, to a renormalized charge density 
in the wire, in the absence of impurity. When the impurity is present, there is in addition a non trivial spreading of 
J> . the charges along the wire. 

The key ingredient in the analysis of |^ is the equivalent of Poisson's equation, which becomes a relation between 
^ . , the electrostatic potential Lp and the charge density: e(/? = u^p. Here, uq is related to the Luttinger liquid constant 
' hy g — {\ + uq/tt Ti vf)~'^^^- The electrostatic potential in turn shifts the band bottom, and thus the total density. 
There follows a relation between the bare injected densities and the true densities : 
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As for the bare densities themselves, they are related with the external voltage sources 
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The hamiltonian including the impurity term reads then, after bosonization 
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where ?; = ^ is the sound velocity. 

To proceed, one defines odd and even combinations of the bosonic field. Only the even field interacts with the 
outside potential, and gives rise to a current. Setting 
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the hamiltonian of interest is 
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where (f>e is a pure right moving field. In these new variables, the boundary conditions (Q) read 
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In going from ffl) to (0), the relation pi?„L 
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, that follows from bosonization, has been 



used. We also have defined pe,o = ^dx<Pe,o- Finally, a mistake in ||] was corrected (see [p^). 

Our goal is to compute the current / flowing through the system as a function of the applied voltage U . In 1^, this 
was accomplished in the case 17=^, using a mapping on free fermions. In this paper, we shall solve the problem for 
general values of g using integrability of the boundary sine-Gordon model Q , . This paper can be considered as a 
sequel - and to some extent a correction - to the work |^,^, where the charging effects were not yet fully understood. 
It is also an extension of the short letter ll6| . 



II. GENERAL FORMALISM 



First, we set e = v = Ti = \ {so vp = g). To treat the interaction term at a; = in an integrable way, one needs to 
chose an appropriate basis for the bulk, massless, right moving excitations, which obey e — p. For g generic, the basic 
excitations can be kinks or antikinks - carrying a pe charge equal to ±1 - and breathers. In the following we shall 
often restrict for technical simplicity to g — j, t an integer. There are then t — 2 breathers. We shall parametrize 
the energy of the excitations with rapidities, ej = rrije^ . Here mj is a parameter with the dimension of a mass; for 
kink and antikink, m± — p, while for breathers, rrij = 2p sin 2(1-1) ' ~ 1, ... ,t — 2. The value of p is of course of 
no importance since the theory is massless, and in the following we shall simply set it equal to unity. The massless 
excitations enjoy factorized scattering in the bulk. At a temperature T, and with a choice of chemical potentials, they 
have densities given by solutions of the thermodynamic Bethe ansatz equations, which we shall generically denote by 
aj (not to be confused with charge densities). 

The key point is that these excitations have also a factorized scattering through the impurity, described by a 
transmission matrix T±±. This matrix depends on the ratio of the energy of incident particles to a characteristic 
energy scale Tb ■ In the following, it is useful to parametrize Tb — e^^ . The modulus square of the transmission 
matrix have very simple expressions; we recall that 

1 + g2(g-iJ(e-eB) 

Finally, we also recall how Tb is related with the bare coupling A H: 
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Ts = (2 sin ^5)^ [Ar(l - g)/2f'^'-'^ . (8) 

To proceed, we start by expressing the boundary conditions in terms of the massless scattering description: 

1 / 

p,{L/2) = — (a+ |r++|V |T+_|^a_ - |T__|^a_ - |r+_|V+ 



^J^Ja+-a^){\T++f-\T+^f)d9. 



(9) 



Here, a± are the densities of kinks and antikinks; one has (7± — nf± where the pseudo energies e± are equal and 
satisfy n = 5^^e±, n = a + a'^ the total density of states of kinks or antikinks (the factor occurs because it is 

the electric charge ^ / dx(t)e associated with the fundamental kinks of the problem). The e's follow from the solution 
of the TBA system of equations 

e,=Tj2N,,^*ln{l + e'^), (10) 

k 

where s{9) = cosh((f-i)e) ' 5 ~ 7' '^^'^ incidence matrix of the following TBA diagram 



1 2 t- 3 

o — o — o — o — O t - 2 



The equations (y^ have to be supplemented by asymptotic conditions sa mje as — s- oo. In (|10[), the chemical 
potential vanishes for all the breathers which have no U{1) charge. For the kinks and antikinks, fi± = T-^j where 
W has to be determined self-consistently (the logic here is that the external potential and the temperature determine 
uniquely the average densities everywhere in the quantum wire. As always in macroscopic statistical mechanics, this 
can be described instead by a distribution with fixed chemical potentials, which is exactly what the TBA allows one 
to handle. By U{1) symmetry, it is known in advance that only the kinks and antikinks have a non vanishing chemical 
potential /i± = ). The filling fractions read then 

The charge density on the left side of the impurity reads simply pe{—L/2) — / (c+ ~ (T-)d9 (it can be simply 
reexpressed in terms of W: Pe{—L/2) = \f\^), so the boundary conditions equation (^) reads 

+ \ |T+ -I') (<^+ - o-)d6 = ^. (12) 

The other key equation in the solution of the problem follows from the charge density drop across the barrier 

V 

Ap = p{x <0)- p{x > 0) ^ g — . (13) 

TT 

Here, p = pn. + PL, and V is the four terminal voltage (that is, the voltage difference measured by weakly cou- 
pled reservoirs on either side of the impurity; it consists of an electrostatic potential drop, plus an electrochemical 
contribution). By following the previous transformations, one finds that Ap = A(ct+ — cr_), and thus dl^ ) reads, 

J\T+^\^{a+~a^)de = g^. (14) 
Finally, the tunneling current / = ^^^^ reads, from (|l^) and ( p^ 
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If T/Tb or U /Tb are large (the high energy, or weak backscattering hmit), the solution of (|T^) is / ((T+ — a-)d9 « 

and thus from / w J{a+ - cr-)d0 = ^. This, once physical units are reinstated, reads I = ^U, the expected 
formula for a spinless quantum wire. 

From the foregoing system of equations, it is now easy to deduce the following identity giving the parameter W in 
terms of the physical voltage and current 

U = 27: (^1-^^ I + W. (16) 

The following relation is also quite useful: 

97r 

V = W-—I. (17) 
9 



III. RESULTS 



The limit 5 — > 1, which describes non interacting electrons, is very simple. In that case indeed, the T matrix 

elements become rapidity independent, and the system of equations can readily be solved to give V = |r_| |^ U, 

I = |T_|_ _(_| Here, the transmission probability is not trivial in general, since, as g —> 1, 9b has to diverge to 
ensure a finite value of the bare coupling A 0,11^ . 

The system of equations determining / can also be solved easily in the "classical limit" g ^ 0, where (this is 
detailed some more in the appendix) 

7^0 ^ sinh(VF/2T) Tb 

^~2.g— — — — -, x = —, (18) 

27r IiW/27vT(^X)I_iW/27vT[^X) AT 

I are the usual Bessel functions, and W follows from (|l6| ). 

Closed form results can also be obtained for g = ^ (see below); besides, except at T = 0, one has to resort to 
a numerical solution of the TEA equations. To tackle the physics of this problem as g varies, we consider first the 
linear conductance at temperature T. In the limit U —> 0, the foregoing system of equations can easily be solved by 
linearization, giving rise to 



1 +1^3^ (t+P 



(19) 



where Go is the linear conductance in the quantum Hall effect problem Q (the numerator of this equation). One of 
the roles of the denominator is to renormalize the conductance from g to unity in the UV region. In the case g = \, 
equation (^9|) can be evaluated in closed form to give 

1 1 IS-Vl/' (1 -I- ) 

^~ 2^1 + ^*'fi + ^)' ^ ^ 

^ ^ 2ttT ^ \2^ 2ttT) 

where ^ is the digamma function. For values g = j,t an integer, G is easily determined numerically by solving the 
system of TEA equations (|l^) , and plotting the soliton pseudo energy back into (^9|) . Curves for various values of g 
are shown in Fig. 1; features entirely similar to those in |^ are observed, although all the curves now converge to the 
same value in the high temperature limit, in contrast with the quantum Hall edges case. The effect of the impurity 
is considerably amplified as q gets smaller, with G getting a discontinuity in the weak back scattering limit as g 0. 
Indeed, letting [/ — > in ([iq), one finds 

C^l,^^_ . . (21) 
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As g ^ 0, G thus becomes a step function, jumping from ^ to as soon as Tb {Tb ~ 2A for g = 0) is turned on, for 
any temperature. 
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FIG. 1. We represent here the conductance as a function of the universal ratio of temperatures T /Tb for several values of 
g = 1/t, t SiU integer. In this domain - which is the easiest to study numerically - G has only a weak dependence on g. These 
curves interpolate between two limiting behaviours: for g = 1, 2nG should become a constant equal to 1/2, while for g = 
2-kG should vanish for any finite value of T /Tb- 



Another simple Umit to study is the case T = 0, where results are far more intriguing. Consider first the classical 
limit: as W is sweeped, one finds that / vanishes while U increases up to ttTb, then goes back to zero, beyond which 
I increases like I = In other words, the system behaves either hke a perfect insulator, or like a perfect conductor! 
This very singular behavior is the 5 — > limit of a multivalued I — U characteristics with regions of negative differential 
conductance [|l6|, that we now study in more details. 

Indeed, the TBA equations can be solved in closed form in the limit T — > . In that limit, (t_ = 0, the integrals 
run only from —00 to a cut off (Fermi) rapidity A, and a+ = n follows from the solution of the integral equation 

n{e)- ^{e-e')n{e')de' (22) 

J-c^ 27r 
while A is determined by the condition that (\{A) = 0, where 

^{9 - e')eX{e')de' ^ 4r - ■ (23) 

oo 2 

In that equation, $ is the derivative of the log of the kink kink S matrix 

sinhTT 



2g-l 
2(l-g) )^ duj 



2cosh^sinh:5;^ 2tt' 

2 2(l-g) 



Since W determines A uniquely (one finds — , where the propagators G are defined below), in what follows 

we will consider instead A as the unknown when T = 0. After a few rearrangements, the relevant equations read now 
(we still set g = \, although t does not have to be an integer here)) 

A /-i_e2(*-i)(e-«B) U 

"(^) . ..r.-nr.-... ^ -, (24) 



and 



l^^2{t-i){e-eB) 27r' 



2(t-i)(e-eB) 

I = I -(^) i + ,2(.-i)(.-..) (25) 



5 



The density n{9) can be computed as a power series in the weak and strong backscattcring Umits, giving rise to 
expansions for the current and the boundary conditions. In the strong backscattcring case one finds ; 



1 = 



while the boundary condition reads 



,^+A-0B^2n(t-l) 



it- 



G+(0) 

In the weak backscattering limit instead, one finds 



^G+(0)" T{n)T{l+n{t-l))^' 



U. 



G+(0)7rt2 



^/^ V(n/t) 



where 



2G+{i) 



tG+{Q) t\t J G+{0) 
Here we have introduced the notations 

G+{lo) = V27rt 



2r(n)r(|+n(i-l)) 

oo 



^A+A-ei5^2n(i-l) 



^ r(n/^) . 4+A-eB\2"(i-i) 



r(n)r(|+n(i-l)) 



i7. 



where 



r(i-i|)r(-ia;5^) 



A = ;^ ln(t - 1) - ^ Ini. 

2 ^ ' 2{t-l) 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



In terms of the auxiliary variable W, the strong and weak backscattering expansions have matching radius of conver- 
gence: either one of them is always converging, and both are at the matching value ^e~^ = 1, where the parameter 



is defined as 



G+(0)- 



Ti, = 2TBe- 



The series can be summed up in the case g — ^ to give 

U -2ttI U + 2-nI 
tan ■ 



(32) 



(33) 
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There is a rich physical behavior hidden in these equations. To investigate it, consider first the behavior of physical 
quantities as a function of W . Curves representing U and — / as a function of W for various t are given in Fig. 2 and 
Fig. 3. 
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FIG. 2. The applied voltage difference U /Tb as a function of the chemical potential difference between solitons and antisoli- 
tons, W/Tb- Observe the remarkable non monotonic behaviour that settles in for small enough values of g. This results in the 
existence of two possible values of W for a given U , and thus in the existence of the loop in the I — U characteristic. 
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FIG. 3. In contrast, the current 7 as a function of W exhibits, once properly rescaled, a very weak dependence on g. 
curves behave asymptotically as W/Tb in the weak backscattering limit. 
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As g 0, the current in the strong backscattering expansion is exactly 0. In the weak backscattering expansion 
meanwhile, it reads 



27r/ 
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hence exhibits a square root singularity at finite value of W (we note that the latter expression can also be obtained 
directly from the result (^8|) by using the uniform asymptotic expansion of Bessel functions for large orders ||l9[] : 



(1 + ^2)1/4^ 



where ry = \/l + + In ■ 



, ). When t is varied, the current evolves from this singular behavior to the simple 

characteristics I — ^ as g ^ 1 (this is easily seen from the integral representations of / and U : and an artifact of the 
variable Tb used throughout, that would have to be rescaled appropriately in that limit to give a non trivial I — U 
relation Q). At fixed g ^ 1, I ^ ^ at large W. 

As 5 ^ 0, [/ in the strong backscattering expansion is simply equal to VF, while in the weak backscattering 
expansion it reads 



U ^ W - {W^ - 7r^Tiy\ 

As 5 — !■ 1 meanwhile, U ~ W. When g varies, U interpolates between these two limiting behaviors, and stops having 
a (local) maximum around t « 4.83. 

The fact that U can decrease as W increases is a direct consequence of the physics in this system. The density 
on the left, pe{—L/2) oc W. An increase in W increases the left density, but it also increases the right density, since 
particles being more energetic, more of them go across the impurity. U is a non trivial function of the densities on 
either side of the impurity, as given by (||). For g large, U behaves essentially as the sum of the densities in ±L/2, 
thus increases when W increases. However, when g 0, U gets dominated by the difference of the densities, and if 
enough particles go across, it can well decrease when W increases. This effect is directly related to the fact that the 
differential conductance Stt^^ does, for g < ^, actually get larger than g for finite values of W an effect first observed 



|] (see Fig. 4). 
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FIG. 4. The rescaled derivative of the current with respect to W at T — 0. Notice the existence of a maximum above the 
weak backscattering limit (equal to 1) for t > 2. This peak of differential conductance becomes more and more marked as 

Consider now / as a function of U: clearly, the existence of a maximum in the curve U (W) will lead to an S-shaped 
I{U). More precisely, consider first the case 3 « 0. Suppose we increase W starting from 0. According to Fig. 2,U 
first increases up to ttTb, then decreases back to zero. W being still finite, / vanishes identically, since it has an 
overall factor of g. Going now to the regime where W becomes infinite, U w W, and ^ ~ ^ ~ the system has 
switched from being a perfect insulator to being a perfect conductor! This is easy to understand in more physical 
terms: as g 0, the kinetic term dominates the Lagrangian, and one might expect that the impurity is essentially 
invisible. However, as 5 — > 0, there is the possibility that a charge density wave might form, getting pinned down by 
an infinitesimal potential, and leading to a perfect insulator jlGt . 

This effect is stable against quantum fluctuations, and for g approximately smaller than g = .2, a "loop" keeps 
being observed in the I—U characteristics. That the current is not a single valued function of U in the region of small 
voltages, leads to the prediction of hysteresis and bistability in the strongly interacting, out of equilibrium regime. 
Although the present calculation is valid only in the scaling regime, this qualitative aspect should survive beyond it. 

The loop is also stable against thermal fluctuations: as is illustrated in Fig. 5 for the case t — 6, it only disappears 
at a finite temperature Tc which depends on g. 
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FIG. 5. We illustrate on this figure the disappearance of the S shape as the temperature is increased. Clearly, the bistability 
is stable against thermal fluctuations in a finite range which depends on g. Here, g — 1/6. 

A semi classical approximation gives Tc = Tb^J'^^^: this formula is not quite correct for values of g < .2, 
but becomes increasingly good as g — > 0. It is quite difficult numerically to determine Tc with a good accuracy: a 
reasonable estimate of this curve is given in Fig. 6. 
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FIG. 6. The "critical" temperature Tc{g) at which bistabihty disappears. Notice the poor quality of the leading semi-classical 
approximation (full curve). 

IV. DUALITY 

For the problem of tunneling between quantum Hall edges, a striking duality between the weak and strong backscat- 
tering limits was uncovered in 1^ at T = 0, and further generalized to any T |20[| . The meaning of this duality was 
that, while the hamiltonian describing the vicinity of the weak backscattering limit is given by (m, the one describing 
the vicinity of the strong backscattering limit can be reduced, as far as the DC current is concerned, to an expression 
identical with (||), up to the replacement of the coupling A by a dual coupling A^, together with the exchange g ^. 
As a result, a duality relation for the current followed 



Here, the dual coupling Ad reads 



2tt 



gl \d,gU, ■ 



(34) 
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[3r(.g)l 







(35) 



The relation (pq) follows from keeping the parameter 



rpll 



Ik 



(36) 



constant]^ while letting .g ^ ^, and using the relation ^ between Tb and the bare coupling in the tunneling hamil- 
tonian. 

For pedagogical purposes, it is probably wise to explain a little more explicitly what the duality means. Consider 
thus a hypothetical current defined non perturbatively by the expression 
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(37) 



It obeys the following duality relation 



1 1 

x' g 



g^i{x,g)- 



(38) 



^In [Q, the duality relation was initially written at constant T'g. While the identities in |^ are algebraically correct, it is really 
Tg that has to be kept constant, since the applied voltage is not left invariant in the duality transformation. 
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Suppose now we did not know the non perturbative expression, but had only access to the small x expansion 



and the large x one 



1 °° / 2 \ ' 



(39) 



(40) 



The duality ( |38| ) could then be deduced from the expansions by say starting from the small x one, setting x — -^,g — ^ , 
and comparing the new expression with the large x expansion. What was done in |^ was to find a similar duality only 
based on the weak and strong backscattering expansions (a non perturbative expression for the current was found 
much later [pT]). 

It is interesting to examine what does remain of this duality in the present case. The IR hamiltonian will behave 
similarly to the case of tunneling between quantum Hall edges, since it is entirely determined by the large A behavior, 
and has no relation with the way the voltage is taken into account. This means that the parameter Tg still has to be 
kept constant in whatever duality symmetry one is looking for. 

There is a quick way to proceed assuming from |^ the relation (p^, which becomes here 



qW f 1 
I(\W,g)^'^-gl[Xa,9W,- 
2tt \ g 

Using this, together with the relation (|l6|), one finds the additional relation 

1~ 



U{\,W,g) = U{Xd,gW, 
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From this it follows that 



I{\U,g) 



U 
2^ 
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(41) 



(42) 



(43) 



For completeness, we can also give a direct proof of this relation. It is convenient first to put the equations in a 
more compact form, namely 



A. 



l-(i-l)^a„AM*-i) 

oo 

z,(i,u,) = A,^a„Af(*-i) 



(44) 



for the strong backscattering limit, and 



A„ 



,2"(i-l) 



t - pn^^ 



2n(i-l) 



for the weak backscattering limit. In (MJ),(E5[) 



_^^_2)n+i2"(*-i) 



T{nt) 



r(n)r(f +n(i-l)) 



T{n/t) 



2 ( ^) r(n)r(i + .(i-i)) 



2 tt/ 



2 U 



t 2T'^ 



(45) 



(46) 
(47) 
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To match with our previous notations, A = G+(i,t)^; however, A in the foregoing equations is determined by the 

^ B 

external voltage, and no reference to or Tg are necessary in its definition. 
It follows from (||) and (H) that 



1 



t 



is {t, u) = -A [ -,tu.^ ] - —i^ ( 7' 



1 



1 



t 



(48) 



where the parameter A is the same in both is and iw Of course, the current is an analytical function of the applied 
voltage, independent of whether one considers the weak or strong backscattering expansions, so the labels s,w can 
actually be suppressed from the equations. It follows that, going back to physical variables, 
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I {X,U,g) ^ —W { -,U - gl { \d,U,- 
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(49) 



Now, W in turn can be expressed in terms of U,I, using the relation (|16[), reproducing ([43[). 

The relation between the current and the applied voltage is implicit in the foregoing equations. It can, however, be 
made explicit by elimination of the parameter A, and we quote here the lowest orders for completeness. In the weak 
backscattering limit one has 



and in the strong backscattering limit 



{t-l){t-2)(3l 



r/32 



(ty)4(i-l) + l + ... 



2(t-l) + l 



[it-l)i2t-l)al + a2] + 

Meanwhile, the parameter A can also be expanded, say in the weak backscattering limit: 

A = /3i(tti)2(*-i)+i + 



(50) 
(51) 



(52) 
(53) 



(54) 



One can directly check the duality relation ( ^ ) on these formulas. Notice that despite the more complex physics, 
which now involves screening, the exponents of the weak and strong backscattering expansions are the same than in 
the fractional quantum Hall case. 

Finally, the duality was extended to finite temperatures in |^^, meaning that formula (^4|) holds at finite 
temperature. Since (16) is still true too, the formula (ES) extends to finite temperature as well. 



V. CONCLUSIONS 



This paper hopefully solves the tunneling problem with a proper treatment of the coupling to the reservoirs, hence 
completing and correcting |^,|l^. We have only treated here the spinless case, but the method extends straight- 
forwardly to the spinfuU case, at least when the spin isotropy is not broken, and the problem maps onto a super 
symmetric boundary sine-Gordon model 

The duality we observed does raise interesting physical questions, in particular concerning the nature of the 
"charges" that tunnel in the weak backscattering limit. We hope to get back to this issue with computations of 
the DC shot noise. 

Acknowledgments: We thank R. Egger and H. Grabert for an earlier collaboration on part of this material, for 
communicating the results of jl2j before publication, and for many illuminating discussions. This work was supported 
by the DOE and the NSF (under the NYI program). 
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APPENDIX A: SEMI-CLASSICAL COMPUTATIONS 



In studying the classical limit, one usually concentrates on the behavior of ej for j finite while g — + 0, that is t — + cxi 
|p3| . This is not sufficient in the study of transport properties, where the knowledge of e±, that is pseudoenergies for 
nodes at the other end of the diagram, are required. The necessary analysis is a bit more complicated then. First, 
it is convenient to introduce the new quantity Yj{9) = e'^^'^^^^'^, and to recast the TBA system, using the identity 



^ + W^) +'^{(^- wh)) - into 



Y, (e + 



2{t-l] 



YAe- 



2{t-l) 



(Al) 



In the limit where g — > 0, we introduce new variables s = j, a = 2^ , and e ^ = l^- , and expand the left and right 
hand sides of equation (Al) to obtain the Liouville equation p3 



The general solution of this equation that is relevant here is ^ 

1 



(2zsin np) 



Jo 



e 2 



(a+is)-ln(2T) 



J- 



62 



{a-is)-ln{2T) 



(A2) 



(A3) 



where Jp are the usual Bessel functions, p — The freedom in the arguments of the Bessel functions a + is - 

X{a — Qfo + i{s — So)) has been resolved by matching with the asymptotic boundary conditions Cj « 2 sin (^ 2(t-i) ) 
as 6* ^ (X). As for the index of the Bessel functions, it is obtained by matching against the result at low energies: 



^ej(-oo)/T 



sinh(j + l)W/2tT 



sinh W/2tT 

We can now compute et-2 by setting s = 1 in the solution (|A3|): one finds 



JpJ-p I 



|^-gfa-ln(2T)j 



It follows that 



The current on the other hand reads 



In the limit i — > oo, this becomes then 



2T 



de 



2T 



(A4) 



d , 1 + e 
m 



-W/2Tg-e/T 



1 + e-2(t-i)(e-eB) d0 1 + e^/^^e-^/'^ 



I = —2smh{W/2T)e 

271 



-e±{eB)/T 



(A5) 



Replacing e± by his classical expression reproduces then the result (18). 



^This is of the general form of solution e"'^ = ^^"laAai^'^'' equation ddx = ^e'^, where A = ^(e^), B = j^(e^). 
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APPENDIX B: LOW TEMPERATURE EXPANSION. 



The remarkable relation 



TT^T^t (Pi 

IiW,T) = IiW,T = 0) + _^3^(w-,r = 0) 



(Bl) 



3 dW^ 

was initially discovered, following a Keldysh expansion of the left and right hand sides, in the context of dissipative 
quantum mechanics in |2^. In (Bl), W is the chemical potential defined in the text - it would coincide with the Hall 
voltage V in the context of the fractional quantum Hall effect ||] . 











TBA 
Weiss 












T/T 


=0.2 











FIG. 7. The dashed line is the order correction to the non equilibrium current as estimated by the equation (Bl). The 
dotted line is the same correction calculated from the TBA at T = 0.2 (it is difficult, for technical reasons, to go below this 
value with enough accuracy). The two curves are in good qualitative agreement: notice that both of them are below the axis 
in the weak backscattering limit. On this figure, t — 7. 



We shall now prove that the current obtained from the TBA does satisfy this relation indeed: as (Bl) involves out 
of equilibrium quantities and the temperature, it provides a very non trivial verification that a Landauer Buttiker 
type approach can safely be applied to integrable quasi particles. 

To start, we recall the general expression for the current (Eq) 



1^1- 
2-K 



1 



1 



1 + ^{^-WjlMT 

where e itself is a function of T . Recall also the value 

e(e = -oo,T)/T _ ^"^^(^ 



1 



l)W/2tT 



sinh W/2tT 



(B2) 



(B3) 



We will only be interested in the terms of order T and in the current: we can therefore drop exponentially small 
contributions, which makes matters considerably simpler. For instance, only the first term in (B2) contributes, and 
the value of e(— oo, T) coincides at this order with its value for T — Q, e(— oo, 0) = Cmin = ■ 

To proceed, we consider the first term in (B2) and assume first that e{0) takes its T = value. The finite T 
corrections (we denote them by (J/^^-*) then entirely arise from a simple generalization of Sommerfeld's expansion in 
the case of free electrons. We use here the same notations as in the appendix of p^. Introducing the function 



H{e) 



1 



r Tb 



2(t-l) 



(B4) 



we find 



1 



27r 1 + e(<=" 



1 °° 1 

-Y- 

27r ^ " 



-W/2)/T 



de'H{e') 



e=W/2 



W 

T 



1 



de 1 + e(^-W2)/T 



(B5) 



13 



Since we neglect exponentially small terms, we can neglect the filling fraction in the first prefactor, and replace the 
bound of integration in the integral by — oo, using the fact that emin < "t- It follows similarly that only terms with 
n even contribute to the series, and therefore, to leading order, 

f-W/2 rp2 



27r 



1 r^'^ d 



(B6) 



The first term is nothing but I{W, T = 0). As for the second, ai is the standard constant of the Sommerfeld expansion 

£2 d 1 ^2 



ai 



2! 



'-de=^- 

del + e^ 6 



At the order we are working, we finally obtain 



/ de 




^ dH 




\de 




'de 


0=A 



(B7) 



(B8) 



where A is the Fermi momentum introduced in the text. One has on the other hand 

dH t-1 



de 



2cosh^(i- \){A-eB) 



To proceed, we must also take into account the changes of e with temperature in the initial expression of the current. 
The leading order correction turns out to be of order then: this gives a second contribution i^/'^-' to the change of 
the current, and shows that there are no crossed terms to this order. 

Neglecting the exponentially small terms as before, the TBA equations for e do not need the introduction of other 
pseudo energies and read 



{e) = e'-T r ^e - 9') In (l + e-(^(«')-M^/2)/T^ 



(B9) 



Integrations by part and Sommerfeld expansion give, as in the study of a leading correction going as T^. We 
can thus write 6(6*, T) = e{0, T = 0) + T^Se, where we find 



Seie) 



9')Se{0')d0' = ~aiT^ 



^{0 - A) 



(BIO) 



This equation is solved by introducing the operator L of |2^]. Calling the integral operator on the left of (BIO) I — K 
(where / is the identity), one has / + L = j-^- Using that ^/ + • = L, it follows that 

5e{e) ^ -aiT"^ 1^'^ ) L{e,A) (Bll) 

Using the value 



}=A 



de 
dd 



W 



determined from |22f| , we find therefore 

e(6»,r) =e(0,T = O)-T2 
Of course, the operator L can be made explicit: 

L(e,e')^L{e' ,e)^(^{e-e')^ 



?,W^2/t 



L{e,A) 



e")(j){e" - e')de" + ... 



(B12) 



(B13) 



(B14) 



14 



The quantity e we use here is related with another quantity e'^ introduced in the main text , and studied in great 



details in by e = — e'ji. In the latter reference, the following identity is established: 



Li9,A) 



2tW 



Using this and integrating by parts, we find 



si(^\t) = t^ 



TT t 



d0L{0,A) 



dH 



So collecting all terms, 



/(T) = /(T = 0) + r^ — W- 



TT ft dH ,.dH 



(B15) 



(B16) 



(B17) 



To conclude, we now turn to derivatives of the current with respect to W at vanishing temperature. The current is 
usually written as 



/(r = 0) 



p{e)H{e)de 



(B18) 



where the density p is given by p = —^-jj-- Using integration by parts, one has 

dW ~ 27r ./^ dW dO 



Taking another derivative, using ( pT^ ) and ( pTsI ), one finds 



1 / 1 



d^I 



dW^ 2ttW V 2t 



dH 



dO 



and thus, comparing with (B17) 



I{W, T) = I{W, T = 0) + t 



TT^T^ d^I 

3 dW^ 



{W, T = 0) 



(B19) 



(B20) 



(B21) 



(this, up to exponentially small terms and higher order analytical terms), thus proving the identity. 

As commented in the main text and in , the differential conductance for g < negative for large enough W/Tb 
(this result does not rely on the Bethe a nsat z, and is a simple consequence of the non linear I — W curve present 
in the Luttinger liquid). It follows from (Bl) that for such values of the current in the fractional quantum Hall 
problem diminishes when T is increased from T = 0, provided W/Tb is large enough. This is a rather counterintuitive 
phenomenon: a priori, one expects that, the larger T , the more energy there is, and therefore the less important the 
backscattering should be. Of course, the current depends on more complex details than the overall energy, and it is 
well possible that W, T, and the non trivial interactions produce an overall less efficient population of quasi-particles, 
even though T is increased. Notice that the current can also decrease when T is turned on now at fixed U, as is clear 
on figure 5. 



To conclude, observe that, using (Bl) together with the duality relation at T = 0, the same relation is found to 
hold to order T^, in agreement with the fact that the duality relation should actually hold at any temperature. 
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